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1. Introduction 

It has been observed by various authors [31 lU El El El EE2 that the Riemann 
problem for certain equations from nonlinear elasticity and gasdynamics cannot 
be solved for all combinations of piecewise constant initial states with shock 
waves, rarefaction waves and contact discontinuities only. For that reason, 
the notion of a delta shock wave and a singular shock wave was introduced 
and employed by authors quoted above, and it was shown that a large class 
of Riemann problems can be solved globally with these additional building 
blocks. The aim of this paper is to study the interaction of one type of these 
new solutions, the delta shock waves, with the classical types of solutions. 

We continue the investigation of the model equation 



(1) {el} 

(2) {e2} 



Ut + (u 2 /2) x = 
v t + ((it - l)v) x = 



initiated in [2]. This system is derived from a simplified model of magneto- 
hydrodynamics. In 0, the authors found a solution for every Riemann problem 
with the initial data (uq,vq) on the left- and (ui, vi) on the right-hand side from 
zero in the following way. 

The eigenvalues of the above system are Xi(u,v) = u — 1, \2(u, v) = u, 
and the right-hand side eigenvectors are ri(«, v) = (0, 1) T , ^("U, v) = (l,v) T . 
The first characteristic field is linearly degenerate and the second is genuinely 
nonlinear. Thus, there are three types of solution. 

(i) When u\ > uq the solution is a contact discontinuity followed by a rarefaction 
wave, 



u(x, t) 



v(x, t) 



x 
t ' 

Ui, 



X < U(jt 

Uot < X < U\t 

X > U\t 



VQ, 

v 1 exp(u - ut), 
«iexp(f - ui), 



x < (no — l)t 

(no — l)t < x < uot 

uot < x < U\t 

X > U\t. 
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(ii) If u\ < uq < u\ + 2, the solution is given in the form of contact discontinuity 
followed by a shock wave, 



u(x, t) 



where v* = v\ 



v(x, t) = < 
2 — uq — u± 



uo, x < ct 

u%, x > ct 

VO, X < (U - l)t 

f*, (uq — l)t < X < Ct 

V\, x > ct, 



2 + ui - uo 

(hi) If uq > u\ + 2 the solution is given in the form of delta shock wave, 



u(x, t) 



Uq, X < Ct 
Ui, X > ct 



( VQ, X < ct} 
v(x,t) = 1 x>ct j+ao(t)D- + ai (t)D + , 

where D~ and D + are the left- and right-hand side delta functions with the 
support on the line x = ct (see below), c = (no + u\)/2, 

st{c-{ui - 1)) st(c-(u - 1)) 

a (t) = , ai(t) = , 

«o — u i uq — u l 

a{t) := ao(t) + ai(t) is called the strength of the delta shock wave, and 

s := c (vi - v ) - ((tt! - l)vi - (u - l)v ) 
is called the Rankine-Hugoniot deficit (see [2])- 

Our aim is to investigate various possible interactions of a solution in one 
of these forms with a delta shock wave. There are five possibilities for this to 
happen. 

Case 1. delta shock wave interact with an another one 

Case 2. delta shock wave interact with a contact discontinuity followed by a 

shock wave from the left-hand side 
Case 3. delta shock wave interact with a contact discontinuity followed by a 

shock wave from the right-hand side 
Case 4. delta shock wave interact with a contact discontinuity followed by a 

rarefaction wave from the left-hand side 
Case 5. delta shock wave interact with a contact discontinuity followed by a 

rarefaction wave from the right-hand side 

We shall always assume that the shock wave or the rarefaction wave starts from 
(0, 0) and the delta shock wave from another point to left or right from zero. 
The initial data are determined by triplets {uq,u\,U2) and (vq,v\,V2). 

We shall now briefly describe what we mean by a solution in the form of a 
delta shock wav e. 

Suppose R\ is divided into finitely disjoint open sets £7j ^ 0, i = l,...,n 
with piecewise smooth boundary curves Tj, i = 1, ...,m, that is Jlj n Oj = 0, 
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UI=i = where fij denotes the closure of fij. Let C(Oj) be the space of 
bounded and continuous real- valued functions on f2j, equipped with the 
norm. Let M(£li), be the space of measures on Oj. 
We consider the spaces 

n n 

C r = l[C(n t ), M r = llM(n i ). 

1=1 8=1 

The product of an element G = (Gi, G n ) G Cr and Z) = (Di, D n ) G .Mr is 
defined as an element D • G = (D\G\, D n G n ) G A4r, where each component 
is defined as the usual product of a continuous function and a measure. 

Every measure on fij can be viewed as a measure on with support in 
This way we obtain a mapping 

m : .Mr -^.M(IJ) 

m(D) = A. + D 2 + ... + Z>„. 

A typical example is obtained when Ml is divided into two regions fix, O2 by 

a piecewise smooth curve a; = 7(i). The delta function <5(a; — 7(t)) G 

along the line x = j(t) can be split in a non unique way into a left-hand side 

D~ G A^(Oi) and the right-hand component D + G A^(Og) such that 

5(x - 7(f)) = a (t)D~ + ai(t).D + 

= m(a (t)Z?" + ai(i)D + ) 

with «o(i) + = 1. The solution concept which allows to incorporate such 

two sided delta functions as well as shock waves is modeled along the lines of 

the classical weak solution concept and proceeds as follows: 

Step 1: Perform all nonlinear operations of functions in the space Cr- 

Step 2: Perform multiplications with measures in the space .Mr- 

Step 3: Map the space .Mr into A4(K+) by means of the map m and embed it 

into the space of distributions. 

Step 4: Perform the differentiation in the sense of distributions and require that 
the equation is satisfied in this sense. 

Note that in the case of absence of a measure part (Step 2), this is the 
precisely the concept of a weak solution to equations in divergence form. 

Following the reasoning in 5 , delta shocks are required to satisfy the condi- 
tion of overcompressibility, meaning that all characteristic curves run into the 
delta shock curve from both sides. It may happen that at a certain point on a 
delta shock curve, overcompressibility is lost. In this case we replace the delta 
shock by a new type of solution which we call a delta contact discontinuity. 
This new concept is introduced in Lemma 1 and Definition 1 below. 

At interaction points our solutions are computed by continuation as contin- 
uous functions of time with values in the space of distributions (as solutions to 
a new initial value problem at the time of interaction) . 

The result of our investigation is that the interaction of a delta shock wave 
with any of three types of solutions to the Riemann problem, (i)-(iii) above, 
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can be described by means of delta shocks and delta contact discontinuities. 
This is summarized in the following theorem. 

Theorem 1. The initial value problem for system (0 with three constant 
states, one of which produces a delta shock, has a global weak solution consisting 
of a combination of rarefaction waves, shock waves, contact discontinuities, 
delta shock waves and delta contact discontinuities. 

The remainder of the paper is devoted to proving this result by going throug 
all possible cases of interaction. 



Case 1. Here uq > u\ + 2, u% > U2 + 2. The speeds of the delta shock waves 
are c\ = {uq + u\)/2 and c\ = {u\ + u 2 )/2. At the interaction point (xo,to), the 
new initial data are 



where 7 denotes a sum of the strengths of incoming delta shock waves. 

Let u = G, v = H + (ao(t)D~ + ai(t)D + ), where G and H are step functions 



and D = ao(t)D + oii(t)D + is a split delta function supported by the line 
x = x + (t - t )c. 
From it follows 



i.e. c = (uq + U2)/2. Since c\ > c > C2, the wave will be the overcompressive 
one, because of uq — 1 > c > U2- Equation © gives 



This equation gives (with a(t) = ao(t) + ai(t))) the following ODE 

a'(t) = c(v 2 - v ) - ((u 2 - l)v 2 - Oo - l)i> ) =: s 6 i, a(t ) = 7. 

The unique solution is given by a(t) = s(t — to) + 7. Substitution of a in the 
equation gives 

a (t) + ai(t) = a(t) = s(t - t ) + 7 

(u - 1 - cs(t - t ) - 7)a (t) + (n 2 - 1 - cs(t - t ) - j)a\ (t) = 

which has a unique solution ao(t), a±(t), since uo ^ u 2 - 

Thus, the result of the first type of interaction is a single delta shock wave. 



2. Interactions with shock waves 





-c[G} + -[G 2 }=0 



c[H}6 + [(G - 1)H]6 + (a' (t) + a[(t))5 

c{a (t) + ai(t))5' + ((uo - l)a {t) + (u 2 - l)ai (*))<*' = 
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Case 2. Suppose that the delta shock wave is given by 
u(x, t) = G(x + a 2 — c\t), 

v(x,t) = H(x + a 2 - at) + P (t)5~(x + a 2 - at) + (3i(t)5 + (x + a 2 - at), 

g_\u Q , x + a 2 < cit ^_ | d , x + a 2 <cit 
\ux, x + ^ycit' \vi, x + a 2 > c\t, 

and that a contact discontinuity coupled with a shock wave is given by 

Vl, x < (ui — l)t 



u 



Ul, X < C 2 t 
Ul, X > C2t 



V*, (Ul — l)t < X < C2t 
V 2 , X > C 2 t, 



where c 2 = (ui + u 2 )/2 < ci (since u% > u 2 , uq > u\ + 2), and v* = ^(2 + ui — 
u 2 )/{2 + u 2 - ui). 

Let us denote by (to, ^0) t ne point where delta shock wave meets the contact 
discontinuity, i.e. this point is the intersection of the lines x + a 2 = cit and 
x = (ui — l)t. 

In the area bounded by the lines x = (ui — l)t and x = uit, the value of u is 
the constant ui. This implies that the delta shock wave runs through it with 
the same speed ci = (ito + u\)/2 as before. Only the values of (3o(t) and Pi(t) 
are changed into, say, Po(t) and Pi(t) due to the existing difference in vo and 
v t . The new strength of the delta shock wave is now si(t — to) + 70, where 

si := ci(v* - v ) - (u - l)(v* - v ), 

and 70 is the strength of the previous delta shock wave in the point (to,xo)- 
Obviously, the new delta shock wave is an overcompressive wave, since uq — 1 > 

Cl > Ul. 

Let us denote by (ti,xi) the point where the new delta shock meets the 
existing shock wave i.e. the point (ii, xi) is the intersection of the lines x + c 2 = 
cit and x = C2t. Let 71 be the strength of the delta shock wave at this point. 
Therefore, we obtain the new initial data 




Vq, X < Xl 
V 2 , X > Xl 



A solution for the new initial data problem will be a delta shock wave with the 
speed c = (uo + «2)/2 < ci, and c is obtained directly from Q in the usual 
way. Again this speed ensures that the obtained wave is an overcompressive 
one, since uq — 1 > c > U2- 

Substituting u and v into (|2j) gives the strength 

a (t) +ai(t) = si(t - h) +71. 

Using this equation and 

(u - 1 - s(t - h) - 7i)a (i) + (u 2 - 1 - s(t - ti) - 7i)di(t) = 0, 

where s := c(v 2 — vo) — {{u2 — l)t>2 — (^o — l)fo), one can find unique ao(t) and 
ai{t), and this proves the above statement. 
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Case 3. Now, uq > U\ > + 2, and the speed of the shock wave c\ = 
(«o + u i)/2 is greater than the speed of the delta shock wave C2 = (ui + u<i)j2. 
Let (to,xo) be the interaction point of these two waves, and let 70 be the 
strength of the delta shock wave at this point. Initial data are now 



where the value of v* is defined as before. 

Similarly to the previous case, the result of the interaction is a single overcom- 
pressive delta shock wave with the speed c = (110+U2) /2 (since uq — 1 > c> 112). 
As before, c is obtained from © and from © one can find ao(t) and a\(t) in 
the same way as above. 

All the way through the contact discontinuity, the delta shock wave has the 
same speed, only ao(t) and a\{t) are changing. 



One can easily see that interaction of a rarefaction and delta shock wave are 
much more complicated. Now we shall deal with this problem. 

As one could see before, the new initial data include a delta function as a 
part. If the right-hand side of u is greater or equal to the left-hand one plus 
2, the new initial value problem can be solved in a simple way as above and 
the result is a single overcompressive delta shock wave. But, when this is not 
a case, the types of admissible solution known so far are not enough to obtain 
a solution. The definition of a new type of admissible solution, called delta 
contact discontinuity, is given below. Its existence is justified by two facts. 
First, a contact discontinuity emerges in the case when one of the characteristic 
fields is linearly degenerate. Second, if a linear equation has a delta function 
as initial data, it propagates along the characteristic lines. These two facts 
inspired the following lemma and the definition of this new type of elementary 
waves. 

Lemma 1. Let the initial data for system (GK^ given by 



where c = {uq + u\)/2 weakly solves the Riemann problem for hllty) . 

Proof. For every ip £ C^°, supp^n {(x, t) : x = (uq — l)t, t > 0} = 0, it holds 




3. Interactions with rarefaction waves 




where uq > u\, but uq < u\ + 2. Then, the function 

( ( Vq, X < (uq - l)t 

tin, X < Ct , , , , 




that 




(v t ,if) + (((u - l)v) x ,tp) = 0. 
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Our aim is to show that this still holds true when it is allowed that suppi^ 
intersects the supports of D~ and D + , i.e. the line x = (uq — l)t. Let us note 
that the condition uq < u% + 2 means that (uq + u\)/2 > uq — 1 so the line 
x = (uq — l)t is on the left-hand side of the shock line x = (no + u\)t/2. 

Equation (^Q) does not contain v, so it is still satisfied. From © we have that 

Vt + ((«0 - l)«)a: = -(UQ ~ 1)(«* ~ ^o)<5 ~ 7^' 

+ (no - l)(v* - v )<5 + (u - 1)7(5' = 
near the line x = (no — l)t. □ 

Usefulness of this lemma will be clear after the interaction of a delta shock 
and rarefaction wave is treated. Then one could roughly see how a solution looks 
like, since the rare faction wave could be approximated with a large number of 
small amplitude non-physical shock waves (see , for example) . 

Another possible use could be in a sort of a wave front tracking algorithm, 
where systems in question posses a solution containing a delta function. 

Definition 1. Consider a region R where u is continuous function and a curve 
r in R of slope Ai(n, v). A distribution (u, v) £ C(R) x T)'{R) is a delta contact 
discontinuity, if v is a sum of a locally integrable function on R and a delta 
function on T which weakly solves (II 121) on R. 

Let us note that the overcompressiveness condition obviously need not hold 
in this case. But, as we already have mentioned before the lemma, a linearly 
degenerate field resembles a linear equation where this type of a solutions exists. 



Also, we shell try to show admissibility of the delta contact discontinuity 
using the entropy and entropy- flux functions for system (1 1121) . Entropy and 
appropriate entropy-flux functions are given by 

r](u,v) =f(u)+g(ve- u )e u 

q(u,v) =e u g(ve- u )u - e u g(ve- u ) + [ uf'(u)du 



= (u - l)r)(u,v) + f(u). 

Substituting the functions u and v in a neighbourhood of the delta contact 
discontinuity support x = (uo — l)t by piecewise constant functions 

' ' x < (no — l)t — e 

(uq — l)t — e < x < (uq — l)t 
(uq — l)t < x < (uq — l)i + e ' 
^v*, (no — l)t + e < x 

where 7 = lim e ^ £ ( v ie + v 2e)- one gets 

T)(u,v) t + q(u,v) x w 0. 
That is, convex entropy condition is satisfied for each entropy function pair. 



u = no, v 



Vis, 
V2e, 



Now, we are returning to the last two cases which covers the rest of possible 
delta shock wave interactions. 
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Case 4. Suppose that a delta shock wave starts from the point —a 2 , a > 0, 
with the speed c\ = (uq + u\)/2 and meets a contact discontinuity followed by 
the rarefaction wave centered at zero. Denote by (io,^o) the meeting point, 
i.e. it is the intersection of the lines x + a 2 = c\t and x = (u± — l)t. As we 
have already seen, the delta shock wave goes through the contact discontinuity 
without speed change (but its strength is changed) and meets the rarefaction 
wave at some point (xo,to)> 

2a 2 2u\a 
t = , x = . 

Uq — U\ Uq — U± 

Let 70 be the strength of the delta shock wave at this point. In order to see 
what could happen, let us approximate the rarefaction wave with a set of non- 
physical shock waves, supported by the lines x = (ui + nr])t, r/ « 1, n £ N. 
(see Fig. 1.) 





x=(u l —l)t 


(«o. v J 


0.04 
0.03 


(/»(f),v(f)/ 

/ / 

X =Zf, / / 


K, v o) 






0.02 


// 

j 1 x=u 2 t 


2 

— a 


v,) 




0.01 


.// 

1 1 ("2. V 2) 
If 


-0.04 


-0.03 -0.02 


-0.01 




0.01 



Fig. 1. 

At least in the beginning, until (u± + rjn) + 2 < uq, the result of successive 
interactions of the delta shock wave with the non-physical shock waves are delta 
shock waves with increasing speeds, with values (uq,vq) on the left-hand side 
and the values on the right-hand side are the values of the rarefaction wave. 
This guide us to look for a curve To := (c(t),t), such that a delta function 
lives on it, c(io) = ^o- The value of u on the left-hand side of T is uq, and 
c(t)/t = x/t on the right-hand side. Inserting the above data for such a curve 
into one gets the following ordinary differential equation 

(3) -c\t)^-n ) + 1 -{(^-) 2 -u 2 )=0, c(t ) = x , 
which has the unique solution 

(4) c(t) = uot — ay^uo — U\)t, t > to- 

Denote by v(t) the value of v\r in the rarefaction wave, v(t) = V2 exp(c(t)/i — 
U2). Substituting expected delta shock wave given by u = G, v = H + 
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ao(t)D To + oti(t)D^ , where G and H are the step functions with discontinuity 
line To, gives 

- c\t)(v(t) - v*)6 + c(t)(a (t) + ai(t))'S 

- c'(t)(a (t) + ai (t))5' + ((^r) v(t) - (n - l)v^j 5 



+ (n - l)ao(t) + 



c(t) 
t 



lj ai(t)\ 5' = 0. 



Since the following ordinary differential equation 

c'(t)(v(t) - v.) - (n - IK + gg ~ lHg . 
c(t) 

has a unique solution (obtained in a simple manner by an integration), the 
strength of the delta shock wave, a(t), is determined. 

Equating the coefficient of 5' with zero, we can compute the two summands 
ao and a± of a. Since 



, eu/2(n - ui) c(t) aJ2(u - «i) 

c{t) = u ° 2Wt — > — = U0 7t — ' 

the obtained delta shock wave satisties rhe right-hand overcompressibility con- 
dition. Overcompressibility condition for the left-hand side is 



(5) {equ4} 



no — 1 > c'(t). 



Now, we have the following two cases. 

(i) If U2 < no — 2, relation (JS} is satisfied trough all the rarefaction wave and the 
resulting solution is a single delta shock wave with the speed c = (no + u-i)j1 
starting from the point (x, t) which is the intersection of the curve To and the 
line x = u%t. 







0.05 






\x={ui— \)t 
\r 




0.04 
0.03 


x=c(t) 

- / / 


t,x) 

/ 

?» 

x = u 2 t <(u — 2)7 








(= \'/ i 

. / / 
// 




K,V„) 




0.02 








\ 0.01 


.// 


(u 2 ,v 2 ) 


2 

—a 






// 




r 






V. 
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-0.01 


0.01 


0.02 0.03 



Fig. 2. 
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After the time t, the solution in this case is given by 



u 



t>t 



x < xt 
x > xt 



\t>t 



V2, X > Xt 



+ 7(£i)<%,t) 



(ii) Suppose that U2 > uq — 2. Then the delta shock wave supported by To is an 
overcompressive wave only until some point (x s ,t s ) lying inside the rarefaction 
wave. (See Fig. 3.) 



\ x= 


K,-l)« 0.05 




(t..x.) 


/ 


V 


0.04 




) \ 

y 


/ 

/ 

/ 
/ t=u 2 t 


("o,"o 


0.03 
X = H, t 

y »•« 


:/ 


/ 

/ 

/ 

/ 




2 

—a 


\ 0.01 

(«i,v,) \ 


:// 


/ 


(u 2 ,v 2 ) 


V 










-0.04 


-0.02 




0.02 


0.04 0.0G 






Fig 


3. 





So, the admissible solution cannot be prolonged along the same curve To- 
Assuming that the rarefaction wave is approximated by a set of small non- 
physical shock waves, the present problem is described in Lemma ^ the right- 
hand side equals uo — 2 + rj, < 77 <S 1, while the left-hand one equals uq. 

In this lemma, the problem is solved by using the new type of a solution - 
delta contact discontinuity. This is exactly what we shall try. That is, suppose 
that the solution consists of the delta function supported by a line Ti : (x — 
x s ) = (uo — l)(t — t s ) going through an area where u has a constant value uo, 
and a shock wave supported by a curve T2 x = C2(i), where C2(t s ) = x s , 
with the left-hand side values no of the function u and the right-hand side ones 
C2(t)/t (a part of the rarefaction wave). All that means that 02(f) should satisfy 
the same equation © as c(t) with the initial data C2(t s ) = x s = c(t s ), i.e. the 
new shock wave is supported by the continuation of the curve Tq. 
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Fig. 4. 

while c^(t) > c(t)/t, the obtained shock 



Since no > c^(t) and «o — 1 < ^(i) wnne c 
wave, supported by the curve T2 is admissible. 

The Rankine-Hugoniot conditions for T2 after the time t = t s imply 

(6) {eS}-c' 2 (t)(v(t) - w*(t)) + - u(t) - (n - l)w*(t) = 0, 

where w* denotes the left-hand 
© simply determines 

u\ ^ t + B 
w*(t) — 



side value of v along the curve T2- Equation 

_ y/i + B ^ eyL p^ UQ -f. _ 2B\ft — u\), 
yt — B 

where B := ay / 2(uo — u\)/2 = ^Jt~ s . 

The value of v between Ti and T2, denoted by w(x,t) has to satisfy the 
equation 

(7) {e6} „ , .. ... , .. 

of the form w(x,t) = V(y), y = x — (uq — l)t. More 
tial data one gets 



w t + (n - l)w x = 0, ion = w*(t). 

The solution to (JJJ) is of the form w(x, 
precisely, using the initial data one gets 

(8) {e7) 

The curve Ti is given by 

(9) {e8} X - (uq - l)t = X s - (U ~ l)t s = (u s - Uq - l)t s = ~t s . 

Substitution of © into © yields V\r 2 = 00 and V(y) G R for (t,x) lying 
between T± and T2, since B 2 = t s . But w G L\ oc C V. 



I exp( (u B + n \J B 2 + y-2B)(B + \J B 2 + y)- Ul ). 

B 2 + yJ 



12 



MARKO NEDELJKOV, MICHAEL OBERGUGGENBERGER 



In order to verify that it is a solution we note that 

\ v*, x < (u - l)t 1 

I W[X, t), X > [Uo — 1)1, X < C2[t) I 

when 7 s ^r 2 is t ne delta function with the strength 7 S obtaining from the initial 
data at (t s ,x s ). Since v(x, t) is constant along the lines parallel to x = (no — l)t 
in a region where u = no it is clear that it is a solution of (J2J). 

In order to see what is going on after the interaction of the delta shock and 
rarefaction wave, one has to consider three different possibilities. 

(a) u <u 2 - 

Then the delta contact discontinuity and shock wave supported by T\ lies 
inside the rarefaction wave since Ti n {(t, x) : x = U2t} = and T 2 n {(i, x) : 
x = u<2,t\ = (actually, c^it) has the line x = u^t as an asymptote, as t — > oo 
(see EI)). 

(b) Uq > U 2 > U - 1. 

Then the delta contact discontinuity stays inside the rarefaction wave and 
the shock wave supported by I^ intersects the line x = U2t at some point (t, x). 

Now, at the point (t, x) we have the new Cauchy problem for (11121) with the 
initial data (uq, w[x, i)), (^2,^2)- Since uq > U2, but uq < 112 + 2, the solution 
is given by 

x < (u + u 2 )(t-t)/2 
x > (u + u 2 )(t - t)/2 

x — x < (no — l)(t — t) 
{u - l)(t - t) < x - x < (u + u 2 ){t - t)/2 
x - x > (u Q + u 2 )(t - t)/2, 

where = t>2(2 + uq — U2)/(2 + U2 — uq). Let us remark that the function w 
equals a constant value along lines with the slope (uo — 1). Denote by T4 the 
shock line x — x = (no + U2)(t — t)/2. This line is a tangent to the curve T2 at 
(x,t). 

Let r% be the line of slope no — 1 starting at (t, x). Since t is a solution to 




not — 2B\ft = U2t, 



we have 



4B 2 _ Au 2 B 



2 



x 



and 



(u -n 2 ) 2 ' ' (n -n 2 ) 2 

Vt + B _ 2 + n -n 2 
w r 3 = —f~ — -v(t) = — v 2 = v*, 

y/t-B 2 + U2-U0 

after the substitution of t and the ending value of the rarefaction wave, v(t) = 
v 2 . So, the function w(t,x) is continuously prolonged by into the area be- 
tween the lines x — x = (no — l)(t — t) (the contact discontinuity line) and 
x — x = (no — U2)(t — t)/2 (the shock curve). 
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The slope of ^ is the same as the one of Ti . That is, there are no interactions, 
and this case is finished. 

(c) u<i < uq — 1. In this case both of Ti and T2 intersects the line x = C2t. But 
as T2 reaches this line at the time t = 4B 2 /{uq — U2) 2 , before the time when T\ 
would intersect it, analysis is the same as in the case (b) (see Fig. 5) below. 





(u (t ,w(x,t)) \ 




(«o.v.) 


Q.l 






(w - v «) 


0.08 




0.06 








\ 0.04 

K. v o) "\ 

\ 0.02 






y 







Fig. 5. 

Case 5. Suppose that a delta shock wave starts from the point (0, a 2 ), a > 
and meets a coupled pair of contact discontinuity and rarefaction wave at some 
point (xo,to). This is possible if uq < u\, U\ > U2 + 2. Suppose that the 
rarefaction wave is centered (starts from (0,0)). 

The point (xo,io) can ^ e easu y found by solving the equations 

u\ + n 2 



x — a 



to 



2 

2a 2 



'-t, x = u\t, i.e. 
2a 2 u\ 

x 



ui — 112 U\ — U2 

In the beginning of the interaction of the rarefaction and the delta shock 
wave the situation is quite similar to the one in the previous case. The solution 
is given by a delta shock wave supported by Tq = {(t,c(t)) : t > to}, where 
c(t) is a solution to 

(10) {e9} _ c '(t)(£|l- U2 )+I((£|l) 2 -n 2 )=0, c(t )=X , 



i.e. 



c(t) = u 2 t + ay / 2(ui - u 2 )i, t > t . 
Equation (|10[) is in fact Rankine-Hugoniot condition for 

The left- and right-hand side coefficients of the new delta shock wave, cto(t) 
and ai(t), can be found in the same way as in the previous case. If (u(t),v(t)) 
is the value of the rarefaction wave, then on the left-hand side of To the new 
delta shock wave takes value (u(t),v(t))\r Q = {c(t)/t,v\ exp(c(i) / 't — u\)) and on 
the right-hand side it equals (^2,^2). Only the overcompressibility condition is 
still in question. The first condition for overcompressibility on the right-hand 
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side is always satisfied, since c'(t) = 112 + a\j2(u\ — Vuj/V&t > U2- For the 
overcompressibility it is necessary that also characteristic lines run into the 
shock from the left-hand side. 

C (t) = U 2 + — ^ < U{c(t),t) - 1 = U 2 + -jz 1, 

i.e. 

a v / 2(ui - u 2 ) > 
2yfl 

This is true until the time t = t s , where 

a 2 (u\— U2) 2/ \f U< 2 \ 

t s = , x s = c{t s ) = a {u 1 -U2)y— + 

Thus, u s = x s /t s = U2 + 2. 
The first case: 112 + 2 > no- 

Then the termination of overcompressibility takes place within the rarefac- 
tion fan and again we are in a position to use the intuition behind Lemma ^ 
i.e. to look for a solution consisting of a delta contact discontinuity supported 
by a curve Y\ and a shock wave supported by some other curve T%. T\ should 
be below T2. 

Ti is the characteristic line of the equation 

v t + (u- l)v x = 

passing trough (x s ,t s ). Using the fact that u = c\(t)/t on Ti = {(t, ci(t)), t > 
t s }, one can find such a function c\ by solving the initial value problem 

c\(t) = ^--l, Cl (t s )=x s . 
The unique solution to the above problem can be easily found 

c 1 (t)=t(-logt + log( Q2 ^ 2 - U2) )+^ + 2 )- 

Using |J2) and the Rankine-Hugoniot condition, the curve T2 = {(c2(t), t) t > t s } 
is uniquely determined by a solution to 

-4W(* - cj f) + - (^) ! ) = 0. *w - 

i.e. 

C 2 (i) = + 0^2(^1 - U2)t. 

One can see that it equals to the function c(t) from the previous case. 
One has to prove that T\ is actually strictly below the curve T2- 
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Since ci(t s ) = C2{t s ) and c^(t s ) = c 2 (t s ) = u a — 1, it is enough to compare 
c[(t) and c' 2 (t), for t > t s . 

i , ,'a 2 (ui — u 2 ) 

c' 1 (t) = u 2 + l-logt + \og' 

= u 2 + 1 + log 



2 

a 2 (ui - u 2 ) 
2t 



, ,., . dy/Ul ~ U 2 

c 2 (t) = u 2 + - 



2t 



c' 2 (t) > c[(t) if 



;il) {elO} V ±= - 1 > log ( V ^ 2> ) , / 



But the last relation is true; one can check it by changing the variables, and 
noticing that ^fy — 1 > log(y), for y G (0, 1). 

Denote by A the region between Ti and T 2 for t > t s . 

The value of u is u(x,t) = x/t inside A. Therefore, Q is satisfied. Now, Y\ 
is the support of the delta contact discontinuity, and we are trying to find the 
value of v in this area. 

First, let v*{t) denote the value of v on the left-hand side of T 2 . The value 
of u there is given by 



1 x c 2 (t) aJ%ui - u 2 ) 
u\r 2 = - t \r 2 = — =u 2 + ^ . 

The values of u and v on the right-hand side of T 2 are u 2 and v 2 , respectively. 
The Rankine-Hugoniot condition for (j2j) gives 



-c' 2 (t)( V2 - V*(t)) + ((U 2 - 1)V 2 - - l) «•(<)) 

Solving the above equation, one gets 



0. 



y/i + B ~ u 2 ) r - 

«*(*) = ^f- t _ B v ^ B = 2 = 

Denote by w(x,t) the value of u inside D. Then w is the solution to the 
linear partial differential equation 



The solution of the above equation is a constant along the characteristic curves 

dx x 

7 := — — = 1, where 7lro is known. 

' dt t ' " 2 

In particular, w tends to infinity near Ti, but in locally integrable fashion 
because v*(t) = 0(1/ (yt — y/t^)) as t —> t s . 

Now, we shall look for an exit of the delta contact discontinuity and the shock 
wave supported by T 2 trough the rarefaction wave. 
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The line x = u^t and the curve x = c\{t) always has an interaction point, 
say (t, x), for uq < u 2 + 2. The equation 

u t = (u 2 + 2)t + t \og{B 2 /t) = \og{t s /t). 

This equation has a unique solution t > t s . 

The next question is whether the curve x = C2(t) intersects the line x = uot 
or not. An intersection takes place, if the equation 

u^t = U2t + 2By/i, i.e. uq — u% = 2B/Vt 



has a solution t > t s . If no < U2, there is no solution. If uq > 112, then 
the solution t = AB 2 /{uq — U2) 2 = 4i s /(uo — U2) 2 is bigger than t s , because 
2 > Uq - u 2 . 

In both cases we have to solve initial data problem for 1)1121) . given by 



u 



Uq, X < X 
Uq, X > X, 



U%. 5 X X 

w(x,x/uq), X > X 



where w is the right-hand side of v in the region A, constant along the charac- 
teristics of 

'x S\ „ dx x 



w t +[- t -l)w x = 0, 7 :_ = 7 -l. 

long the line x = uot the slope of these characteristic curves is uq — 1. Thus we 
may continue the solution to the left of x = uot as a delta contact discontinuity 

f u„ a;-s<(tio-l)(t-t)l , . 

«(*,,) = »o, KM) = I to(M) _ a _ . > K _ 1)(t _ (1 j+7A- fc( »-ix«,. 

where w(t,x) is again constant along the lines with slope uq — 1. Denote by 
T3 the line x — 5 = (uo — l)(t — t). There is no further intersection with the 
original contact discontinuity along the (parallel) line x = (uq — l)t. In the case 
u o < v>2i the solution is complete (See Fig. 6). 
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Fig. 6. 
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In the case uq > U2 we still have to consider the region above the intersection 
point (x,t) of x = C2(t) with x = u$t. In this case 

u 2 < u < u 2 + 2 

and we can connect a constant left-hand state Uo to the constant right-hand 
state U2 by by a shock wave in u. 

This shock wave supported by T$ has speed (110+112) /2 and actually is tangent 

to the line x = C2(t) at the intersection point (x,t). It follow a classical contact 

discontinuity starting from the point (x, t) with speed uo — 1, supported by the 
line T3 and this connects the region when v = w(£, t) has been determined by 
the initial data along the line x = u^t. The value of v between I^ and T5 is 
v* = (2 + u - u 2 )/ (2 + U2 - uo). (See Fig. 7). 
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Fig. 7. 



The second case is no > u s = U2 + 2. Then there is no bifurcation of the delta 
contact discontinuity supported by To- After To intersects the line x = uot at 
(ti,xi), say, the solution can be continued into the region (uo — l)t < x < uot 
by a simple delta contact discontinuity on the line : x — x\ = (uq — l)(t — 1\), 
where u has the constant value uq and v has the value and v 2 on the left and 
right-hand side, respectively, and a constant strength delta function is placed 
on the line T3. This concludes investigation of all possible cases. 

References 

[1] Bressan, A., Hyperbolic Systems of Conservation Laws, Preprins S.I.S.S.A., Trieste, Italy. 
[2] Colombeau, J. F., Elementary Introduction in New Generalized Functions North Holland, 
1985. 

[3] Hayes, B. T. and Le Floch, P. G., 'Measure solutions to a strictly hyperbolic system of 

conservation laws', Nonlinearity 9, 1547-1563 (1996). 
[4] Keyfitz, B. L., 'Conservation laws, delta shocks and singular shocks', In: M. Grosser et 

al: Nonlinear Theory of Generalized Functions, Research Notes in Math., Champman 

Hall/CRC, 1999. 



18 



MARKO NEDELJKOV, MICHAEL OBERGUGGENBERGER 



[5] Keyfitz, B. L. and Kranzer, H. C, 'Spaces of weighted measures for conservation laws 
with singular shock solutions', J. Diff. Eq. 118,2, 420-451 (1995). 

[6] Korchinski, D. J., Solution of a Riemann Problem for a 2 x 2 System of Conservation 
Laws Possessing No Classical Weak Solution, PhD Thesis, Adelphi University, Garden 
City, New York, 1977. 

[7] P. D. Lax, Hyperbolic Systems of Conservation Laws and the Mathematical Theory of 

Shock Waves, SIAM, Philadelphia, 1973. 
[8] Nedeljkov, M., 'Delta and singular delta locus for one dimensional systems of conservation 

laws', Math. Meth. Appl. Set. 27, 931-955 (2004). 
[9] Nedeljkov, M., 'Second delta locus - interaction of delta shock with shock waves', Preprint. 
[10] Oberguggenberger, M., Multiplication of Distributions and Applications to Partial Dif- 
ferential Equations, Pitman Res. Not. Math. 259, Longman Sci. Techn., Essex, 1992. 
[11] Oberguggenberger, M. and Wang, Y-C, 'Generalized solutions to conservation laws', 

Zeitschr. Anal. Anw. 13, 7-18 (1994). 
[12] Tan, D., Zhang, T. and Zheng, Y., 'Delta-shock waves as limits of vanishing viscosity for 
hyperbolic systems of conservation laws', J. Diff. Eq. 112, 1-32 (1994). 



